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The New South Wales population of the Mountain Pygmy-possum Burramys parvus was estimated 
at 500 adults in 8 km 2 of boulder-heath habitat during surveys in the I980’s. In 1989, this estimate 
was increased to 1312 adults based on further surveys in four of the identified habitat patches. 
However, further research indicated that population density varied greatly between habitat 
patches and that the original extent of habitat and the revised population estimate had been 
overestimated. Because of the high degree of uncertainty regarding the total population size of 
B. parvus in New South Wales and the relative value of each habitat patch, including those in ski 
resorts, it was necessary to re-evaluate these estimates. Results based on data collected between 
1996 and 2001 indicated an adult population of around 6I3±92 in 1.85 km 2 of boulderfield 
habitat.This population was distributed in scattered colonies within a modelled bioclimatic range 
of 444 km 2 of suitable climate, based on modelling of the then known southern Kosciuszko B. 
parvus population. However the notional extent of suitable climate has been greatly enlarged by 
the recent discovery of a northern Kosciuszko population. Establishing which boulderfields to 
target for future surveys within the enlarged potential climate envelope is not simple. No strong 
association with any habitat characteristics were identified by habitat analysis, however, deeper 
boulderfields with a variety of boulder sizes and reasonable vegetation cover were generally 
preferred. Changes in habitat characteristics across the elevational gradient partly explained the 
difficulty in modelling habitat preference. 
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Introduction 

The Mountain Pygmy^possum Burramys parvus is a small 
(35g), terrestrial burramyid marsupial. Described from 
fossil remains (Broom 1896) it was first discovered as an 
extant species when found in a ski lodge in Victoria in 
1966 (Anon 1966a, b) and first caught in New South 
Wales in Kosciuszko National Park in 1970 (Calaby et al 
1971). Surveys in the early 1970’s and 1980’s indicated 
populations were largely confined to patches of boulden 
dominated habitat and associated shrubby heathland in 
the Australian Alps from Mt Bogong to Mt Hotham in 
Victoria and in southern Kosciuszko National Park in New 
South Wales (Calaby et al 1971, Gullan and Norris 1984, 
Caughley 1986, Mansergh et al 1989). 

Caughley (1986) sampled a range of heath types in 
Kosciuszko National Park above 1500 m elevation with 
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and without boulders. The majority of B. parvus was caught 
in rocky areas with a diverse association of shrubby heath 
that included Mountain Plum-pine Podocarpus lawrencei. 
The total population size of B. parvus was estimated at 
around 2000 individuals, which included some 500 adults. 
Caughley’s estimate was based on an average density of 
2.6 ±1.7 B. parvus per hectare, with an adult to juvenile 
and male to female ratio of around 1:3 (i.e., adult density 
0.65/ha), in approximately 8 km 2 of habitat. This estimate 
was derived from 10 survey sites each sampled for 3 nights 
between February and April during a full moon which, 
Caughley argued, sampled the majority of the trappable 
population. At the time of Caughley’s surveys it was not 
known that B. parvus adults begin to hibernate as early 
as late February and most individuals are hibernating by 
late April. 
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Broome (1992) obtained average densities (between 
November and March) of two to three times Caughley’s 
estimates in four of the areas previously used to derive 
estimates. Using a comparison of the most reliable data 
set (from Mt Blue Cow 1986 - 1989), Caughley’s estimate 
was multiplied by a factor of 2.27 to produce a revised 
estimate of 1312 adults in New South Wales (327 males: 
985 females) (Broome and Mansergh 1989). However, 
subsequent trapping by Broome in several additional areas 
surveyed by Caughley (1986) suggested that the four sites 
monitored by Broome (2001a), and one additional site 
on Mt Kosciuszko appeared to have high abundance of B. 
parvus (average density 11.05 adults/ha) compared to many 
of the other sites within the alpine area. The relative value of 
such sites would be underestimated by assigning an average 
density to them, and equally, applying densities derived from 
them to the total 8 km 2 of habitat mapped by Caughley 
(1986) would inflate the total population estimate. 

Much of the habitat mapped by Caughley (1986) was 
identified from high level (1:24,000) aerial photographs. 
Preliminary ground truthing between 19864996 indicated 
that some of these areas were quite unsuitable for B. parvus 
and suggested that the extent of habitat had also likely been 
over-estimated. However, it appeared there may also be 
additional sites, for example, north of the air photo coverage 
near Gungartan, that were not included in Caughley’s 
estimate. Because of the high degree of uncertainty regarding 
the total population size of B. parvus in Kosciuszko National 
Park and the lack of knowledge concerning the relative 
value of each habitat patch, including those within ski 
resorts, it was necessary to re-evaluate these estimates. 
The aim of this study, conducted between 1996-2001, was 
to re-trap and survey habitat areas identified by Caughley 
(1986), identify and survey additional sites, and assess 
habitat characteristics, population abundance and density 
to produce a revised population size estimate and habitat 
map for B. parvus in New South Wales (Broome et al 2005). 
However in 2010-11 three new colonies of B. parvus were 
discovered near Cabramurra, well to the north of the area 
considered in the original study, and outside the predicted 
area of suitable climate based on the southern colonies 
alone (Schulz et al 2012a, Broome, Bates and Shi unpubl. 
data). Climatic predictions were therefore re-run to examine 
the implications of this discovery for the enlarged NSW 
distribution of B. parvus. 

Methods 

Identification of habitat and trapping of 
B* parvus 

Breeding populations inboulderfields (technically periglacial 
blockfields, blockstreams, moraines, or solifluction deposits, 
Rosengren and Peterson 1989, Galloway 1989, or tors and 
boulder scatters around rocky summits) were targeted for 
trapping between November and mid-December. Burramys 
parvus breeds at one year of age, and unless the first litter is 
lost, usually have one litter per year (Mansergh and Scotts 
1990, Broome 2001b, Heinze et al 2004). Broome (2001b, 


unpubl data) found that the highest numbers of adults 
were trapped early in the breeding season in November- 
December and that populations were concentrated around 
well-defined boulderfields at this time of year. Lactating 
females suckle young in the nest during late December to 
January, and frequent trapping during this period may affect 
the survival of young. Dispersal occurs between February 
and April. Adult males and a high proportion of juveniles, 
especially males, disperse and migrate from breeding sites 
to occupy surrounding areas of heathland where boulder 
patches are often small or dispersed (Mansergh and Scotts 
1990, Broome 2001b, Heinze etal. 2004). These habitats are 
critical for the species, but are often difficult to delineate. 

Boulderfields often occurred in clusters, and although 
a range of boulderfield types within each cluster were 
sampled, what was anticipated to be the best habitat in 
each cluster (based on previous field experience) was 
always trapped to maximise the number of individuals 
detected. As a population estimate was a primary aim 
of this study, larger boulderfields had a correspondingly 
higher number of traps. However, density of trapping 
effort was relatively constant. Broome (2001b) found that 
B. parvus individuals frequently move distances of 500 m 
or more between boulderfields within a habitat cluster 
on a nightly basis. Therefore boulderfields less than 500 
m apart that were on the same aspect were trapped and 
treated in the analysis as a single habitat unit. Trapping 
commenced in December 1996. Following helicopter 
reconnaissance of the sites identified by Caughley (1986) 
in June 1997 the bulk of remote sites were trapped 
between November and December 1997. Additional sites 
identified from aerial photos, or on the ground during 
trapping and habitat surveys were trapped between 1998 
and 2001. For the four sites that were monitored each 
year (Broome 2001a, unpubl. data), a year between 1996 
- 2001 was randomly chosen to provide data for this study. 

Trapping was conducted for 3 or 4 nights at most sites. Traps 
at sites trapped for 3 nights were baited with walnuts. Traps 
at sites trapped for 4 nights were baited with chocolate 1 
on the first night to enable collection of scats for dietary 
analysis that were unadulterated by bait. Walnuts were used 
on nights 2-4 to provide a preferred food source. Burramys 
parvus are easily trapped. Caughley (1986) conducted 
trapping over 3 nights during her surveys and found similar 
numbers of possums were trapped each night and no 
evidence of trap shyness or differential trappability between 
age or sex of animals. Trapping on the four sites monitored 
by Broome (2001a) was conducted over 3 nights between 
1987 - 1994 and 4 nights from 1995 - 2001. Chocolate was 
not as attractive as walnut as a bait and fewer animals were 
generally caught on the first night when chocolate was used 
(Broome unpubl. data), however, results generally indicated 
that few additional animals were trapped on the third night 
of baiting with walnuts. Most survey sites trapped in the 
current study, apart from the long-term monitoring sites, 
had insufficient numbers of animals captured to estimate 
population size from capture-recapture models. However, 


' The attractant of the chocolate bait appeared to be the smell and it was seldom ingested. When collecting scats we now use half a cotton make-up pad 
soaked with walnut oil as an attractant for the first night of trapping. 
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modelling on the four monitoring sites by Broome (2001a, 
unpubl.data) over the 15 years from 1987-2001 (using the 
heterogeneity model in program CAPTURE, White et al 
1978, 1982) usually indicated population closure (81% 
of 87 models), high capture probabilities and on average 
an 11% difference between the population size estimates 
(N-hat) and numbers trapped (N). Similarly the average 
difference between the population estimate ± 2 standard 
errors was 15%. In the current study, abundance and 
densities are derived from the numbers trapped, and a final 
adjustment is made to estimate population size by adding 
11 %, the average difference between N and N-hat on the 
monitoring sites. 

Mapping of habitat and measurement of 
habitat variables 

All boulderfield patches were delineated on aerial photos 
enlarged to a scale of 1:2000 or 1:5000. These patches 
were ground-truthed during the trapping and habitat 
assessment surveys. Distortion in the aerial photographs 
was corrected by using reference co-ordinate points marked 
on the photography or by comparing the photographs with 
geometrically corrected satellite imagery. Digitised and 
geometrically corrected boulderfield polygons were stored 
as an Arc View Shapefile (Lourens et al 2002). Eighty-four 
survey sites, including 186 of 605 digitised boulderfield 
polygons, were trapped and surveyed for habitat variables. 

Elevation and aspect were determined for each site from 
a grid based digital elevation model (DEM) and aspect 
surface previously developed for Kosciuszko National 
Park. Percent snow cover was calculated as the percentage 
of each site covered by snow on each of 8 Landsat TM 
images taken between July and November 1996 and July 
and September 1998 (Lourens et al 2002). 

Measurement of a range of habitat variables at each of the 
trapping sites commenced in February 1998. Each site was 
assessed at a site level and at multiple plots within each 
site. Site level variables included slope, the quantity of 
visible Bogong Moths Agrotis infusa (ranked on a scale of 
0-3 from absent to abundant), the distance to the nearest 
watercourse, whether watercourses within the boulderfield 
were absent, intermittent or permanent, continuity with 
other boulderfields (scaled 1-4 from completely isolated 
in grassland to well connected by other boulderfields and 
heath), a visual assessment of percentage vegetation cover, 
the percentage occurrence of 9 vegetation communities 
surrounding the boulderfield (to a distance of 20 m), a list of 
vegetation species present and a visually assessed complexity 
rating. Complexity was initially assigned a scale of 1-3 (little 
boulder structure and low vegetation cover and diversity to 
highly structured, deep boulders and good shrub cover) but 
this was changed to a scale of 1-5 during data preparation 
because many of the values given in the field fell between 
the categories (1-2 or 2-3). Plot level variables included 
percent cover of vegetation, substrate (% cover of rock and 
soil), shrub height (Podocarpus and other shrubs), % cover 
of eight key food plant species (Smith and Broome 1992) 
and grass, presence of all other plants (species for shrubs, 
categories for herbs e.g., Ranunculus spp., daisies, ferns etc.), 
the number of boulder layers (0=rock protruding from soil, 
1=one layer - rock sitting on soil, 2=two layers - substantial 


rocks sitting on other rocks, 3 = three layers - more than 

2 layers of rock but ground still visible and 4 - more than 

3 layers, ground not visible), boulder depth (the vertical 
height from the ground surface to the average height of the 
top layer of boulders, if depth was indeterminable (generally 
greater than 1.5 -2 m) a maximum assessable value was 
determined and 50 cm was added for data analysis), ground 
visibility (scored 0 - no ground visible to 5 - ground only, 
no rock) and boulder size - the number of boulders in each 
of 5 size categories (<20 cm, 20-50 cm, 50-100 cm, 1-2 m, 
>2 m) that could be seen within or touching the plot. Plots 
were circular with a lm radius. The number of plots was 
proportional to the size of the boulderfield, ranging from 10 
on the smaller boulderfields to 90 on the large boulderfields. 
Plots were evenly spaced on the boulderfields following 
contours and extended 10 m beyond the edge to sample the 
surrounding habitat. 

Climatic site profiles and distribution 
modelling 

Thirty-five bioclimatic parameters were generated 
for each trapping site using the BIOCLIM module of 
ANUCLIM (Nix 1986). These were used in habitat 
regression models (below). The Kosciuszko 25 metre 
DEM was used to model five climatic variables (maximum 
temperature, minimum temperature, rainfall, radiation 
with rainfall and evaporation) using the twelve standard 
Australian climatic surfaces provided with ANUCLIM. 
A second climatic profile was produced for all sites 
where populations of B. parvus were recorded during 
this study. The BIOMAP module of ANUCLIM was 
used to produce a predictive map based on modelling 
of six variables; annual mean temperature, maximum 
temperature of warmest period, minimum temperature of 
coldest period, annual precipitation, precipitation of the 
wettest period. The full range (0-100) was considered the 
most appropriate measure as some of the largest known 
populations of B. parvus inhabit the highest peaks in 
Australia (Mt Kosciuszko and Mt Townsend), suggesting 
that the potential climatic distribution of the species 
extends well above the altitudinal availability of habitats. 

Data analysis 

Field collected habitat characteristics were measured both 
at the site level, and at multiple plots within each site. 
Those collected at multiple plots were reduced to single 
variables for each site, or, for characteristics with more than 
one measurement category, to summary statistics based 
around levels of measurement (e.g. percentage occurrence 
of three different boulder size classes across plots). Methods 
of reduction for each variable are listed in Table 1. Further 
reduction in the number of analysed variables was achieved 
via principal components analysis for boulder layers, boulder 
counts, vegetation cover, complexity and BIOCLIM 
attributes. The “complexity” variable was a loosely defined 
visual assessment of the diversity of vegetation and boulder 
structure present within a site. It was therefore not included 
in regression analysis. Various methods were attempted 
to reduce complexity of floristic data from the presence/ 
absence of species at each site. However, the generally 
sparse and low diversity vegetation of boulderfields meant 
that no method identified any distinct patterns. For 
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Table I . Methods used to reduce plot-based variables to site-based variables for analysis. 


Variable 

Initial measurement 

Reduction method 

% Vegetation 

% at n plots 

Classes: % plots <20%, 20-80%, >80% 

Rock 

% at n plots 

Average for n plots 

Soil 

Not analysed: reciprocal of rock 


Height P. lawrencei 

cm at n plots 

Average for n plots 

Height other shrubs 

cm at n plots 

Average for n plots 

Boulder depth (cm) 

cm at n plots 

Average for n plots & %> 1 m 

Visibility 

Not analysed 


Boulder layers 

Number of layers at n plots 

Classes: % plots 0,1,2,3+ layers 

Boulders: <20cm 

Count at n plots 

Average for n plots 

20 - 50cm 

Count at n plots 

Average for n plots 

50 - 100cm 

Count at n plots 

Average for n plots 

1 - 2m 

Count at n plots 

Average for n plots 

>2m 

Count at n plots 

Average for n plots 

Cover key plants 

(each of 9 species) 

% cover at n plots 

Average for n plots 


example, the best principal component of plant species 
occurrence explained only 19% of the variation, and only 
two explained more than 5%. A similar lack of clear pattern 
was obtained using multi- dimensional scaling. Plant species 
occurrence data was therefore represented in analysis only 
by floristic diversity. Variables included in regression models 
were; area (ha), aspect, average boulder depth, average % 
rock cover, average Podocarpus lawrencei height, average 
shrub height, average % snow cover, bioclim PCI, blocks 
water course, Bogong Moth count, boulder count PCI, 
boulder count PC2, B. parvus presence, abundance and 
density, continuity, % cover of each of eight key food plants, 
grass, and all vegetation combined, elevation, floristic 
diversity, layers PC land PC2, slope, surrounding proportion 
of snowgum woodland, tall heath, short heath, damp bog, 
raised bog, herb/grassfield, valley bog, fen and feldmark, 
vegetation cover PCI, water distance. 

Interactions were expected between many variables; 
altitude and aspect in particular were expected to affect 
vegetation attributes and snow cover. Scatterplot matrices 
were used to identify correlated variables, and regression 
analysis was used to test the strength of relationships. 
Number of traps (and trap nights) were correlated with 
total possum captures (0.52), but not possum density 
(0.08) or detection of possums at a site (0.35). Ordinarily a 
correlation between number of traps and trapping success 
could indicate sampling bias or methodological problems. 
However, sites with less traps were smaller boulderfields, and 
could therefore not hold more traps. Intensity of trapping 
(traps/unit area) was not correlated with captures (-0.18). 
Some trapping “bias” is evident due to the monitoring 
grids, which are known to be large populations, and which 
contain some of the highest numbers of traps due to annual 
trapping programs. The correlation between trap success 
and number of traps was significantly lessened by excluding 
the four monitoring sites (0.35), and it is a reasonable 
assumption that the number of possums captured at survey 
sites was representative of the size of those populations at 
the time of the survey. Exploratory regression analyses or 
t-tests were used to identify variables that were most likely to 


influence presence, abundance (individuals/site) or density 
(individuals/ha) of B. parvus . The number and density of 
B. parvus and the more common small mammal species, 
Bush Rats Rattus fuscipes , Dusky Antechinus Antechinus 
swainsonii and Broad-toothed Rats M astacomys /uscus, 
were positively correlated across the sampled boulderfields. 
Other small mammal captures were therefore not included 
in regression analyses, as it was thought that the correlations 
were likely due to similar needs for food resources and 
shelter and any direct interactions between possums and 
these species would have a negative effect on possums, 
rather than increasing their abundance. 

Multiple regression analyses were performed in the 
generalised linear modelling module of Statistica 6 to 
determine what factors influenced; 1) presence of B. 
parvus in a boulderfield, 2) abundance of B. parvus within 
a boulderfield, 3) density of B. parvus within a boulderfield. 
Variables included in models were those identified in 
exploratory analysis as being potentially different between 
sites with and without B. parvus , and those that exhibited 
reasonably strong correlation with the abundance or density 
of B. parvus. Models initially included all such variables, and 
were reduced via step-wise removal of the least significant, 
until only significant variables remained (at the 0.05 level). 
All variables were then individually re-added to the model 
and were included if they were found to be significant. 

Results 

Trapping of B* parvus 

Eighty-four trapping sites were clustered into seventeen 
geographic units based on average distance between 
boulderfields using eastings and northings in JMP 3 (SAS 
Institute). These were increased to 20 clusters on the basis 
of topography by separating Dicky Cooper from Whites 
River, Muellers from Townsend, Etheridge from Kosciuszko, 
and including Byatts Camp in the Townsend cluster rather 
than with Wilkinsons Valley (Figure 1). A compilation of 
trapping results for these clusters are presented in Table 2. 
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Figure I. The distribution of boulderfields in the Kosciuszko Plateau region of Kosciuszko National Park. Numbers refer 
to boulderfield clusters described in Table 2. 
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Re-assessment of Mountain Pygmy-possum in Kosciuszko National Park 


A total of 338 B. parvus were captured at 50 sites. Low 
numbers of B. parvus (five or fewer females) were trapped 
in half of the habitat clusters. The highest numbers were 
obtained at Charlotte Pass, Blue Cow-Guthega, Townsend- 
Byatts and Kosciuszko (Table 2). Numbers of B. parvus 
trapped fluctuate from year to year, but the population 
estimates on the four monitoring sites in 1997, when a large 
proportion of the sites were trapped, were not significantly 
less than the 15 year averages (Broome 2001a, unpubl. 
data). Numbers trapped on the monitoring sites during the 
years randomly selected to include in this study did not 
differ from the means by more than 2 individuals, with the 
exception of number of females at Charlotte Pass, which 
was 5 above average (this would make only a 1% difference 
to the percentage totals in Table 2). 

Other small mammal species loosely associated with B. 
parvus in the boulderfields included M. fuscus (0.25) 
and A. swainsonii (0.30) which were trapped in low 
numbers at approximately half of the clusters. Agile 
antechinus Antechinus agilis was rarely recorded, while 
R. fuscipes (0.51) occurred in 80% of the clusters (Table 2). 
Numbers of these species on the four monitoring sites were 
also as high or higher than the 10 year averages in 1997, 
except for low numbers of R. fuscipes on the Paralyser site. 

Trapping results from sites surveyed by Broome (unpubl 
data) in previous years are presented in Appendix 1, 
results from Caughley’s surveys in Appendix 2 and details 
of the trapping results from this study in Appendix 3. 

Abundance and density of B. parvus were only moderately 
correlated (0.59), making modelling of each necessary. 
No single variable was a good predictor of B. parvus 
presence, abundance or density. Complexity, which was 
not included in modelling analysis (see methods), was 
the only variable obviously correlated with B. parvus 
presence (0.56), abundance (0.48) and density (0.44). 
However, although not directly modelled, significant 
factors in the final regression models (Table 3) reflected 
the importance of characteristics that constituted the 
“complexity” variable. The importance of boulder structure 
was indicated by PCI of boulder layers and PC 2 of boulder 
count being significant in all models. PCI of boulder layers 
distinguished boulderfields that had more than two layers 
(negative values, associated with possums) from those that 
did not (Table 4), while PC2 of boulder count distinguished 
boulderfields dominated by boulders greater than lm in 
diameter (positive values), from those that had a more even 
distribution of boulder sizes (negative values, associated 


Table 3. Significant variables from regression models. 


Variable 

Presence/ 

absence 

Abundance Density 

% Snow cover 

n.s. 

(+) 0.0021 

n.s. 

Floristic diversity 

(+) 0.0318 

(+) 0.0021 

n.s. 

Boulder count PC2 

(-) 0.0402 

(-) 0.0013 

(-) <0.0001 

Layers PC 1 

(-) 0.0024 

(-) 0.0031 

(-) <0.0001 

Layers PC2 

n.s. 

(-) 0.0131 

n.s. 

Boulder count PCI 

(+) 0.0103 

n.s. 

n.s. 

Veg. cover PC 1 

(+) 0.0146 

n.s. 

(+) 0.0062 


Table 4. Weightings for significant principal components. 


Category 

PCI 

PC2 

Eigen value 

2.73 

1.32 

Percent 

54.59 

26.45 

Cum percent 

54.59 

81.03 

BC <20cm 

0.55 

- 0.12 

BC 20-50cm 

0.55 

-0.25 

BC 50-100cm 

0.54 

-0.06 

BC l-2m 

0.32 

0.58 

BC >2m 

0.06 

0.76 


Eigen value 

1.95 

1.32 

Percent 

48.79 

33.2 

Cum percent 

48.79 

81.99 

PL0 

0.48 

-0.49 

PL 1 

0.54 

0.41 

PL 2 

-0.39 

0.63 

PL 3+ 

-0.57 

-0.45 


Eigen value 

1.92 


Percent 

63.84 


Cum percent 

63.84 


PV <20% 

-0.71 


PV 20-80% 

0.41 


PV >80% 

0.57 



BC= Boulder class, PL Proportion of layers, PV Proportion 
veg class. 


with B. parvus ). PCI of boulder count was also significant 
in the presence model, indicating that presence of boulders 
of all sizes was important to the presence of possums. PC2 
of boulder count was significant in the abundance model, 
which is more difficult to explain, as negative values 
(associated with possums) identified sites with more three - 
layered plots, but also more plots with no layers, possibly 
indicating a diversity of habitats within a site. Vegetation 
characteristics contributed to all models. There was a 
significant positive correlation between floristic diversity 
and possum abundance and presence, and a significant 
positive correlation between the degree of vegetation cover 
and density and presence of B. parvus. The percentage snow 
cover was significant in the abundance model, with sites 
with high numbers of possums tending to have high snow 
cover. None of these models was a particularly good fit to 
the data (R 2 = 0.2849 presence, 0.3253 density and 0.3274 
abundance: Figures 2 and 3). However, deeper boulderfields 
with a variety of boulder sizes and reasonable vegetation 
cover are obviously implicated as preferred habitat. 

The expected effect of altitude on other variables was 
readily apparent, and helps in interpreting some of the noise 
associated with the modelling results. Floristic diversity was 
a significant factor in presence and abundance models. 
However, there was a clear negative association between 
altitude and floristic diversity (-0.75), and altitude and 
vegetation cover (-0.60), yet some high altitude sites 
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Figure 2. Model fit for possum presence. Line represents 
whole model prediction of possum presence. Circles 
are the actual proportion of occurrence of possums at 
sites within a particular prediction class (reduced to 10 
categories for ease of presentation. Error bars are +/-1 
standard error 

(Kosciuszko and Townsend in particular) have high 
numbers of possums despite their low vegetation diversity 
and cover. This is probably explained by the fact that plants 
are important food sources, but at higher altitudes where 
plants were less abundant, another important food source, 
Bogong Moths Agrotis infusa, were more abundant than at 
low altitudes (correlation 0.53), and could offset lack of 
other foods (Smith and Broome 1992, Gibson 2007). 

Estimation of population size 

Calculating an unadjusted population estimate (N) by 
multiplying the average density estimates by the total area 
of habitat within each boulderfield cluster (Table 2) will 
almost certainly produce a high estimate of the population 
size in most clusters. This is because the boulderfield 


areas not trapped were of lower complexity, i.e. had few 
layers, low vegetation cover or were shallower than those 
selected for trapping. This was particularly evident for 
areas such as Whites River, Muellers Peak, Kosciuszko, 
Summit Road and Paralyser, where there was a relatively 
large amount of untrapped, low quality boulderfield. 
Adding 11% to the resultant cumulative estimate of 677 
in Table 2 (to adjust for the difference between numbers 
caught (N) and the expected population estimate (N-hat) 
and adding a confidence interval of 15%, see methods) 
produces a very optimistic estimate of 750± 112 adult 
B. parvus in Kosciuszko National Park. An alternative 
method is to apply an average density estimate to the 
untrapped portion of the habitat and add this to the 
total for the trapped sites. Use of the unweighted average 
density over all trapped boulderfields (338 animals / 86.7 
ha = 3.90/ha) results in a high estimate of (338+383) + 
11% = 800+120. However, the median density (2.18/ 
ha) is likely to produce a more realistic estimate because 
less weight is given to the highest densities. This produces 
a total of (338 + 214) + 11% = 613+92 adult B. parvus 
(395 females and 218 males) in southern Kosciuszko 
National Park between 1996-2001. 

Very conservatively our results show that 24% of the 
southern Kosciuszko population was within the boundaries 
of the Charlotte Pass and Perisher-Blue ski resorts, but this 
may have been as high as 39% (Table 2). 

Prediction of habitat distribution and 
modelling of habitat variables 

The modelling of potentially suitable climates based on the 
known colonies of B. parvus surveyed between 1996-2001 
in southern Kosciuszko (South Ramshead to Gungartan 
Pass) predicted areas above 1600 metres in a limited 
region of southern Kosciuszko that included areas of the 
Cascade ranges and around Mt Jagungal, where possums 
have not been recorded (Figure 4). Boulderfields on Mt 




Figure 3. Model fit for possum abundance (a) (individuals/site), (b) possum density (individuals/ha). Solid line is the 
predicted occurrence with 95% confidence limits shown as dotted line. 
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Figure 4. Climatic predictions based on known populations of Burramys parvus in NSW.The “south” prediction is based 
on populations studied in the survey component of this study the “all” prediction includes the recently discovered 
northern population. 
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Jagungal have been trapped twice (in 1987 and 2009) 
with no success, while no welhdefined boulderfields have 
been noted on the Cascades. The total area predicted as 
suitable climate was 44,416 ha. Within that area 185 ha 
of boulderfield was mapped, with 75 ha total area known 
to contain possums. 

In 2010-11 three new colonies of B. parvus were 
discovered at elevations below 1660 m near Cabramurra 
in the northern section of Kosciuszko National Park 
(Schulz et al. 2012 a, Broome et al 2012, Broome, Bates 
and Shi unpubl. data). Inclusion of successfully trapped 
boulderfields near Cabramurra dramatically increased 
the potential distribution of B. parvus in NSW, as the 
new potential climatic envelope (207,864 ha) includes 
a number of previously unsurveyed boulderfields 
as far north as the Bogong Peaks and Brindabella 
Ranges (Figure 4). Boulderfields identified from aerial 
photographs and satellite imagery within this broader 
region will be targeted for future surveys, and may 
contain populations of possums not considered in the 
population estimate presented here. However, the 
density of boulderfields within the wider potential 
distributions is markedly lower and more patchy than 
in the southern distribution. 

Discussion 

We estimated the total number of adult B. parvus in 
Kosciuszko National Park during these surveys (1996- 
2001) to be around 613 ± 92 but optimistically numbers 
may range up to 860. Despite the fact that Caughley’s 
(1986) estimate was based on a lower mean density of 
adult B. parvus than we found, our estimate is closer to 
the 500 adults estimated by Caughley (1986) than to the 
1312 revised estimate (Broome and Mansergh 1989). 
This is because the heterogeneity of boulderfield habitats 
is far greater than previously recognised, and the total 
area of B. parvus habitat, particularly high quality habitat, 
identified by this study was significantly less than that on 
which Caughley based her estimate. It is evident from the 
current survey that the densities of B. parvus in the four 
monitoring sites were higher than the average or median 
densities, hence the inflated estimate of Broome and 
Mansergh (1989). Although our approach was to map and 
survey boulderfields and not include surrounding areas 
of heath with scattered boulders as Caughley did, it is 
evident that the extent of heath and boulderfield habitat 
as mapped by Caughley (1986) was overestimated in some 
areas. Mapping largely from high level (1:24,000) aerial 
photographs, Caughley (1986) included areas that had 
only surface rock and large areas of grassland, particularly 
on the western scarp of Abbott Ridge including Mt 
Townsend, where half the habitat (4.35 km 2 ) was thought 
to occur. This led to the assumption that over half of the 
population occurred in that area (rather than 12% as 
indicated by this study), well away from any development 
activity. Further, it was assumed that the population of B. 
parvus in the proposed Blue Cow ski resort represented 1% 
of the numbers in Kosciuszko National Park (Caughley 
1983, 1984, National Parks and Wildlife Service 1985) 
rather than 14-20% indicated by this study. 


The largest populations of B. parvus occurred in the two 
ski resort lease areas, with the next highest numbers on the 
peak of Mt Townsend and the western fall of Abbott Ridge 
(cluster 7 Townsend-Byatts Camp) and Mt Kosciuszko. 
However, except for the very peaks of these mountains, 
densities were generally quite low in the latter two areas, 
reflecting extensive areas of relatively low quality habitat. 
Many of the boulderfields on the western slopes of Abbott 
Ridge e.g., at Byatts Camp and those in Wilkinsons 
Valley were characterised by very large boulders and only 
one male possum was trapped at those sites. Vegetation 
cover was low and boulder depth quite shallow in 
the extensive boulderfields on the west side of Mt 
Kosciuszko. The habitat analysis indicated that preferred 
habitats were characterised by deeper boulderfields with 
a variety of boulder sizes and reasonable vegetation 
cover. Floristic diversity was important in explaining 
the presence and abundance of B. parvus but was not 
significantly correlated with density. This was because 
some of the highest densities occurred in boulderfields 
on the mountain peaks, including boulderfields on the 
summits of Mt Townsend, Muellers Peak, Kosciuszko and 
Etheridge Ridge, where there was low vegetation diversity 
or cover but good boulderfield structure and high numbers 
of Bogong Moths. Bogong Moths form a large part of the 
diet of B. parvus in summer and attain highest densities 
in boulderfields at high elevations (Mansergh et al 1990, 
Smith and Broome 1992, Broome 2001b, Green 2010). 

Previous research has shown that high quality B. 
parvus habitat is comprised of deep boulderfields at 
high elevations with associated shrubby heathland 
characterised by Podocarpus lawrenceii and Leucopogon 
spp., and high numbers of Bogong Moths (Mansergh 
and Scotts 1990, Mansergh and Broome 1994, Broome 
2001b). At Mt Blue Cow the population uses a mosaic of 
boulderfield patches and shrubby heathland, with some 
males travelling and returning on a nightly basis up to 
1.5 km between boulderfield patches and females up to 1 
km from low to high elevations to access Bogong Moths 
(Broome 2001b). Shrubs provide the principal autumn 
dietary source of seeds and berries (Mansergh et al 1990, 
Smith and Broome 1992, Gibson 2007) and also provide 
habitat for caterpillars, beetles, spiders and millipedes, 
arthropod dietary items of secondary importance for B. 
parvus. Unvegetated boulderfields (screes) were found by 
Green (1988) to have low numbers of invertebrates in 
pitfall and litter samples. 

The Kosciuszko Peak boulderfield was surveyed a number 
of times between 1990-1998 (Appendix 1). In 1990, 77 
individuals were trapped during the combined February 
and March surveys (15 adult females, 12 adult males 
and 50 juveniles). In April 1992, the heaviest individual 
ever trapped, an 82 gram female fattened for winter 
hibernation (Broome and Geiser 1995) was trapped at this 
site. The site is unvegetated, but during these surveys high 
numbers of Bogong Moths were aestivating at the site. 
However, only 2 females were trapped in December 1993, 
a near or above average year for B. parvus populations 
on the monitored sites (Broome 2001a) but few Bogong 
Moths were present. It is likely that the populations of B. 
parvus in this cluster and other sites such as Townsend and 
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Muellers are highly dependent on the migratory Bogong 
Moths and these populations may vary in size according 
to conditions that affect Bogong Moth populations in 
their lowland breeding grounds. Possums may move 
between the rocky peaks and the shrubby slopes below, 
as was found at Mt Blue Cow (Broome 2001b), however 
at Kosciuszko and Muellers nearby shrubby heathland 
is sparse. Habitat areas with good vegetation cover 
and diversity that support relatively large populations 
include those at Summit Road, Charlotte Pass, Spencers 
Creek, Paralyser and Blue Cow. Possums in these areas 
would have a greater potential to diversify their diets so 
population stability is likely to be greater at these sites. 
These shrubby sites were also where R. fuscipes occurred 
in high densities. 

High numbers of B. parvus were not detected by Caughley 
(Appendix 2) on many of what we have found to be key 
sites because either the principal, high quality parts of the 
site were not trapped (generally the high altitude rocky 
peaks), or trapping was conducted after many individuals, 
particularly adults, would have commenced hibernation. 
It was not known at the time of Caughley’s surveys that 
adult B. parvus hibernate and become untrappable from 
late February and juveniles from March or April (Broome 
and Geiser 1995, Walter 1996, Kortner and Geiser 1998, 
Broome 2001b). Mt Stillwell, Muellers and the Paralyser, 
for example were trapped by Caughley in late April and 
the main Charlotte Pass boulderfield in May. Unaware of 
the overwinter hibernation strategy of B. parvus , Caughley 
(1981) reported that the reason no B. parvus were trapped 
at Muellers in late April was not obvious. While Summit 
Road was trapped by Caughley in March, only juveniles 
were captured and the subsites trapped on Mts Blue Cow 
and Townsend included some of the mkbelevation sites, 
but not the peaks where highest densities occur (Broome 
2001b, Appendix 3). 

Our 1996-2001 surveys resulted in a small extension of the 
known range along Kerries Ridge at the northern end of 
the then known distribution area (noted but not mapped 
by Caughley due to a limit in aerial photo coverage). 
However, BIOCLIM modelled a much larger area north of 
the Main Range and Kerries Ridge as suitable climate for 
B. parvus. Broome (unpubl. data) trapped Mt Jagungal in 
1987 and found no evidence of B. parvus or habitat suitable 
for the species on the mountain peak. Boulderfields on the 
western slope of Mt Jagungal were trapped in 2009, again 
with no success. The area was extensively burnt in the 2003 
wild fires and all mature P lawrencei perished. Although 
substantial boulderfields are not evident from small-scale 
aerial photos, ridges in the Brassy, Dargals and Grey Mare 
Ranges would be worthy of further investigation, as would 
the Chimneys Ridge south of the Thredbo River (Figure 
2). Caughley surveyed around Mt Selwyn in October 1980 
with negative results. Evidence of B. parvus presence (cross- 
sectioned P lawrencei seeds, Mansergh and Broome 1994) 
was detected at some small patches of boulderfield that were 
not trapped in this survey viz. Farm Creek between cluster 
19 and 20 (Figure 1) and on the slopes of Disappointment 
Ridge, (south of cluster 2 and north of cluster 20, Figure 1). 
A female was trapped at one of the Disappointment Ridge 
sites in February 1999 (Glen Sanecki pers. com.) and two 


males were trapped by K. Green in February 2003 and April 
2004 at the northern end of Kerries Ridge, approximately 
3 km north of cluster 1. However, it was not expected that 
substantial additional areas of habitat would be found, 
because of the restriction of known breeding populations 
in the southern Kosciuszko region to elevations above ca. 
1600m elevation, variant slightly on latitude and aspect 
(Caughley 1986, Broome 2001a,b, this study), although 
captures of individual animals had been made as low as 
1400 m. As our results have illustrated, any new areas 
can be rapidly assessed by experienced observers using the 
complexity rating. 

Recent discoveries of an additional three B. parvus colonies 
between Mt Jagungal and Cabramurra (Snow Ridge, Rough 
Creek, Boltons Hill - Happy Jacks Valley) in the northern 
Kosciuszko region resulted in the capture of 29 adults (15 
females, 14 males) and 47 juveniles (22 females and 25 
males) during November 2010 - January 2011. Surveys at 
the same sites during October - November 2011 and 2012 
resulted in the capture of 160 adults (124 females and 36 
males) (Schulz et al. 2012a, b, L. Broome, H. Bates, H. Shi 
unpubl. data). Potential boulderfield habitat between Mt 
Jagungal and the Brindabella Ranges has been digitised 
from aerial photo coverage. Surveys were conducted in two 
of these identified sites during October 2011 (Byron and 
Kiandra), at Farm Ridge and Toolong Range in December 
2012 and Bogong Peaks in February 2013 with no captures 
of B. parvus. Trapping in January 2012 on a spoil dump 
deposited by the Snowy Hydro-electric Scheme in the 
1950s near Guthega (Guthega Adit Camp) on the northern 
banks of the Snowy River close to cluster 20 (Figure 1) 
resulted in the capture of an additional 5 adult females 
and 4 adult males (Schulz et al 2012b). At this point, 
predictions using densities from the successful sites and 
mapping of potential habitat suggests that locating more 
than an additional 250 adult females and 80 adult males is 
unlikely (L. Broome, H. Bates, H. Shi unpubl. data). 

Provision of a point-in-time estimate of total population 
size is necessarily approximate because changes occur 
due to the influence of climatic trends, population 
sizes of colonies can fluctuate independently due to 
local conditions (e.g., Broome 2001a) and uncertainty 
occurs in the accuracy of estimates. There has also been 
the additional, unexpected discovery of new colonies. 
However, the total population on the Kosciuszko plateau 
(optimistically about 395 adult females and 218 males 
south of Mt Jagungal estimated from this study in 1996- 
2001, with perhaps an additional 250 adult females and 80 
males in newly discovered colonies north of Mt Jagungal) 
is small. Most local populations of breeding females 
are much smaller than 20, which has been cited as the 
minimum number below which demographic stochasticity 
is a real threat (Soule 1987, Hanski 1991). McCarthy and 
Broome (2000) suggested that because of the strength 
of density dependence, populations of B. parvus greater 
than 15 females are relatively safe from extinction from 
stochastic processes under (then) current conditions. 
However, only a third to a half of the clusters meet this 
criterion, unless there is substantial movement between 
them. Genetic analysis of hairs from each capture will be 
carried out to determine the degree of isolation of clusters. 
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In 1994 the total population size across the species 
range was estimated at about 2600 adults, 1300 in 
Victoria and 1300 (an overestimate as indicated by 
this study) in NSW (Broome and Mansergh 1989, 
Mansergh and Broome 1994). An additional population 
was discovered at Mt Buller in 1995, increasing the 
Victorian estimate by 300 adults (Heinze and Williams 
1998). Heinze et al. (2004) presented a revised estimate 
of 2056 adults (1570 females, 486 males) in which the 
population at Kosciuszko was suggested to be around 
250 adult females. This was due to declines in the 
southern Kosciuszko colonies that began in the late 
1990’s. By 2009, the numbers in the four colonies in 
the southern Kosciuszko distribution that exceeded 25 
females (clusters 7, 10, 17, 20 in Table 2) had declined 
by between 36-77%, with an overall decline of 43% 
and a population estimated at around 355 adults (225 


females and 130 males, Broome et al. 2012). Trapping 
on the southern Kosciuszko long-term monitoring sites 
in November 2011, 2012 and 2013 indicated a recovery 
at Charlotte Pass and Kosciuszko and a partial recovery 
at Mt Blue Cow (currently 50% of 1986-99 numbers). 
Mt Townsend was not trapped, but given unprecedented 
high numbers on the other monitoring sites (Summit 
Road and Paralyser) and at Gungartan, Gungartan Pass 
and Whites River, it is probable this colony has also 
recovered. These recoveries follow a return to four years 
of high rainfall. It is therefore likely that the current 
population on the southern Kosciuszko plateau is similar 
to that estimated in this study. In round figures, with the 
addition of predicted numbers from colonies north of Mt 
Jagungal, the adult population on the Kosciuszko plateau 
in 2011-2013 is likely to be in the order of 950 adults ( ca. 
650 females and 300 males). 
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Cootapatamba w'fallA.B.C 12-14/12/97 0.71 45 3 135 0 0 0.00 

Cootapatamba X,YZ 12-14/12/97 0.74 60 3 180 II 2.71 

Swampy Plain A 4-7/12/00 0.50 25 4 100 2 I 6.04 

Swampy Plain B 4-7/12/00 0.30 35 4 140 3 0 9.97 
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APPENDIX 4 


Re-assessment of Mountain Pygmy-possum in Kosciuszko National Park 



High quality habitat at Mt 
Blue Cow (1700-1984m) 
in 1996: diverse Podocarpus 
lawrencei heath land overlaying 
deep boulders with a high 
component of small to 
medium sized rocks.This site 
was burnt in 2003. Photo, L. 
Broome. 



The peak of Mt Blue Cow 
(this photo) and the peaks of 
Mt Townsend and Kosciuszko 
are unvegetated but can have 
high numbers of Bogong 
Moths, a favoured food item 
for the Mountain Pygmy- 
possum. L and R Broome with 
moth trap. Photo, M Schroder 



High quality habitat at 
Charlotte Pass on the slopes 
below Mt Guthrie (1750 - 
1920m). Photo, L. Broome. 
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H. Shi and W. Zhongming 
measuring rock structure 
among diverse Podocorpus 
lawrencei heath land at the 
Paralyser 1760- 1900m 
(extends to 2054m on Mt 
perisher). Photo, L Broome. 



Mountain Pygmy-possum 
adult female trapped at the 
Paralyser Photo by R. Plum 
and C. Corkran. 



The steep western slopes 
of Abbott Ridge with Mt 
Townsend 2209m (left) and 
Abbott Peak 2159m (centre). 
Photo, I. Pulsford. 
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Re-assessment of Mountain Pygmy-possum in Kosciuszko National Park 



Volunteers trapping on 
the western slopes of Mt 
Townsend. Note the large 
size of the rocks; densities of 
Mountain Pygmy-possums are 
low in these areas. Photo, L. 
Broome. 



The peak 2228m (left) 
and western slopes of Mt 
Kosciuszko. This area can 
support high numbers of 
Mountain Pygmy-possums but 
annual numbers can be quite 
variable. Photo, L. Broome. 



The Whites River boulderfields 
(1670 m) are bisected 
by a road and power line 
easement. Possum numbers 
were low here despite good 
habitat structure. Photo, M. 
Dawson. 
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